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Abstract 



We consider a lattice formulation of the four dimensional N — 1 Wess-Zumino model that 
uses the Ginsparg- Wilson relation. This formulation has an exact supersymmetry on the 
lattice. We show that the corresponding Ward-Takahashi identity is satisfied, both at 
fixed lattice spacing and in the continuum limit. The calculation is performed in lattice 
perturbation theory up to order g 2 in the coupling constant. We also show that this Ward- 
Takahashi identity determines the finite part of the scalar and fermion renormalization 
wave functions which automatically leads to restoration of supersymmetry in the continuum 
limit. In particular, these wave functions coincide in this limit. 
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1 Introduction 



Recently, there have been several attempts to study supersymmetric theories on the lattice 
PQ-[H] (f° r recent reviews and a complete list of references, see jl]). The major obstacle 
in formulating a supersymmetric theory on the lattice arises from the fact that the su- 
persymmetry algebra, which is actually an extension of the Poincare algebra, is explicitly 
broken by the space-time discretization. Without exact lattice supersymmetry one might 
hope to construct non-supersymmetric lattice theories with a supersymmetric continuum 
limit. This is the case of the Wilson fermion approach for the N = 1 supersymmetric 
Yang-Mills theory j3j where the only operator that violates the N = 1 supersymmetry is 
a fermion mass term. By tuning the fermion mass to the supersymmetric limit one re- 
covers supersymmetry in the continuum limit (see Ref. [5] for numerical studies along this 
approach) . 

The strategy of most recent studies is to realize part of the supercharges as an exact 
symmetry on the lattice. This exact supersymmetry is expected to play a key role to restore 
the continuum supersymmetry without (or with less) fine-tuning of the action parameters. 
These ideas apply to theories with extended supersymmetry where the lattice theory is 
realized by an orbifolding construction U\ . Another approach is based on writing the 
theory in terms of twisted fields jHJE]- The connection between twisted fields and Kahler- 
Dirac fermions is emphasized in ^U] and recently in [TT] . 

In this paper we consider the N = 1 four dimensional lattice Wess-Zumino model 
introduced in Refs. [121 C3j and studied in [T3] where it was shown that it is actually 
possible to define a lattice supersymmetry transformation which leaves invariant the full 
action at fixed lattice spacing. This transformation is non-linear in the scalar field. The 
action and the transformation are written in terms of the Ginsparg- Wilson operator and 
reduce to their continuum expression in the naive continuum limit a — > 0. In [H] the algebra 
of this lattice supersymmetry transformation was studied and the closure of the algebra 
was explicitly shown to g 2 order. This is a necessary ingredient to guarantee the request of 
supersymmetry. It was also argued that the existence of this exact symmetry is responsible 
for the restoration of supersymmetry in the continuum limit. In this paper, we derive the 
Ward-Takahashi identity (WTi) associated with this lattice supersymmetry transformation 
and show how in the continuum limit one recovers the WTi associated with the continuum 
supersymmetry transformation. This will be done in lattice perturbation theory up to 
order g 2 . An outcome of this approach is the calculation of the lattice renormalization 
wave function for the scalar and fermion fields. 

The paper is organized as follows. In Sec.|2]we briefly review the N = 1 four dimensional 
lattice Wess-Zumino model based on the Ginsparg- Wilson fermion operator, and show how 
to build up a lattice supersymmetry transformation which is an exact symmetry of the 
lattice action. In Sec. El we derive the WTi and we explicitly check the simplest one, the 
one-point WTi at one-loop. A second and more interesting WTi, relating the boson and 
fermion two-point function, is analyzed at g 2 order in Sec. 0] Here it is shown that this 
identity is exactly satisfied on the lattice. In Sec. E]we verify this WTi in the continuum 
limit and determine the finite part of the lattice renormalization constants which allow to 
identify the continuum invariant theory. Technical details of the fermion Dirac operator 
and the tadpole cancellations are presented in Appendix A and B, respectively. 
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2 The Wess-Zumino model 



We formulate the lattice Wess-Zumino model by introducing a Dirac operator which satisfies 
the Ginsparg- Wilson relation jT3] 

l5 D + D l5 = aD l5 D . (1) 

This relation implies the existence of a continuum symmetry of the fermion action which 
may be regarded as a lattice form of the chiral symmetry ^B] and protects the fermion 
masses from additive renormalization. As shown in Ref . ^3] , using this Dirac operator it is 
possible to introduce a local action which is chiral invariant and where the fermions satisfy 
the Majorana condition. Moreover, in order to keep as much as symmetry as possible, 
the bosonic kinetic operator must be written in terms of D. The lattice action for the 
Wess-Zumino action reads 

Swz = Sq + S int , (2) 



with 



S ° = E - l^i)" 1 ^ - l^DiA + BD\B) 

X * 

2 v 2 ' 2 v 2 ' J 

f 1 1 



(3) 



X 



^9[F(A 2 - B 2 ) + 2G(AB)} } 



(4) 



where A, 5, F and G are real scalar fields and \ is a Majorana fermion which satisfies the 
Majorana condition 

X = X T C (5) 
and C is the charge conjugation matrix which satisfies 

C T = -C , CC ] = 1 . (6) 

Moreover, our conventions are 

Ci.C- 1 = ( 75 ) T • (7) 
The operators D\ and .D2 which enter in Sq are related to the operator D in (0) by 

A = ^Tr(D) , D 2 = ^(7^) . (8) 

Our analysis is valid for all operators that satisfy Eq. ((TJ), however, in the following we will 
use the particularly simple solution given by Neuberger [T7j 



where 



and 



A = ^(v; + v M ) - ^v;v M 



(10) 



-(4>(x + a/t) — 0(x)) 
a 

—((f)(x) — (f)(x — aft)) 



(11) 



are the forward and backward lattice derivatives, respectively. Substituting Eq. (JHJ) in 
Eq. (jHJ) one finds 



D x = - 

a 



:i + yv;v M ) 



A = ^(v*+v^- 



7A' (I 2 ) 



The Ginsparg- Wilson relation JI} implies the following relations for .Dx and D 2 



D\ - D 2 



■A 



(13) 



and 



;i-fi)-^ 



1 n2 

2 



■A 



(14) 



Before concluding this section we list the propagators of the lattice perturbation theory 
for the scalar and fermion fields: 



(AA) 

(FF) 
(AF) 

(xx) 



where 



(BB) = -M- 1 (l--D 1 y l 



(GG) = -M- l D 1 = - 
a 

(BG) = mM' 1 
((l-^A^A+m) 



-M-\l- a - Dl )-'Dl 



-M~\{1 - |A) -1 A - m) 



A* 



-a^-a)-^™ 2 



(15) 



(16) 



and the Ginsparg- Wilson relation (JH)) has been used to rewrite the auxiliary fields propa- 
gators. Despite the appearance of the operator (1 — |Z}i) _1 , there are no would be doublers 
and the propagators are regular (see appendix A). 



2.1 The supersymmetric transformation 

As discussed in jT2], So is invariant under a lattice supersymmetry transformation which 
is obtained from the continuum one by replacing the continuum derivative with the lattice 
derivative D 2il - On the contrary the interaction term Si nt breaks this symmetry because of 
the failure of the Leibniz rule at finite lattice spacing |]Q . In order to discuss the symmetry 
properties of the lattice Wess-Zumino model one possibility is to modify the action by 
adding irrelevant terms which make invariant the full action. Alternatively, one can modify 
the supersymmetry transformation in such a way that the action (J2J) has an exact symmetry 
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for fixed a. In jHj it has been shown that the full action (J2J) is invariant under the following 
supersymmetry transformation 

SA = ex = xe 

SB = -ie>y 5 x = -ixi^ 

S X = -D 2 (A - i lb B)e — (F — i l5 G)e + gRe 

5F = eD 2X 

5G = i eD 2l5X , (17) 

where R is a function depending on the scalar fields and their derivatives that can be 
determined in perturbation theory imposing the invariance of the Wess-Zumino action 
under (|17|). 

By expanding R in powers of g, 

R = rW + g Rj?) + ... (18) 
and imposing the symmetry condition order by order in perturbation theory, one finds 

R W = ((i _ ~D 1 )~ 1 D 2 + m)~ x AL (19) 

with 

AL = -^=^2(AD 2 A — BD 2 B) — D 2 (A 2 — B 2 ) + 2ry 5 (AD 2 B + BD 2 A) — D 2 {AB) } (20) 
and 

R(n) = _ V2((l - ^D 1 )~ 1 D 2 + m)-\A + i^B)R^ , (21) 

for n > 2. Notice that the operator ((1 — ^Di)~ 1 D 2 + m) -1 is precisely the free fermion 
propagator and that the transformation (|17|) . like the function R, is non-linear in the scalar 
fields. Indeed, using flT^j) and (|2*T|) one sees that the expansion (fTHj) can be resummed and 
R is the formal solution of the equation 

[(1 - ^D 1 )~ 1 D 2 + m + V2g(A + i lb B)]R = AL . (22) 

Notice that, in the limit a — > the transformation (|17|) reduces to the continuum super- 
symmetry transformation, since AL vanishes in this limit. Indeed, AL is different from 
zero because of the breaking of the Leibniz rule at finite lattice spacing. 

In ^3] it has been shown that the algebra associated with the lattice supersymmetry 
transformation (fT7|) closes. The existence of this exact symmetry should be responsible for 
the restoration of supersymmetry in the continuum limit. 

In the following sections we will prove that the Ward-Takahashi identity (WTi) derived 
from this lattice supersymmetry is exactly satisfied at finite lattice spacing. We will perform 
a one-loop analysis though the procedure can be generalized to higher loops. We will also 
discuss the a — > limit. 

3 One-point Ward-Takahashi 

The WTi is derived from the generating functional 

Z[$, J]= J T?$exp -(S wz + Sj) (23) 
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where Sj is the source term 

s j = j* - $ = X) \ Ja a + Jb b + Jp f + Jg g + r ■ ( 2 4 ) 

Using the invariance of both the Wess-Zumino action and the measure with respect to the 
lattice supersymmetry transformation ()17|) . the WTi reads 

(J*-^)j = 0, (25) 

with <5$ given in (|17jl . 

We begin with the simplest WTi which is obtained by taking the derivative with respect 
to fj and setting to zero all the sources 

(D 2 (A - i l5 B)) + (F) - z 75 (G) - g(R) = . (26) 

The order 0(g) of this Ward-Takahashi identity is 

(D 2 (A - i^B))^ + - il5 (G)^ - g(Rp-))W = , (27) 

where the notation (O)^ indicates the n— order (in g) contribution to the expectation 
value of O. From Eq. ([T5|) it is easy to see that all the terms of the WTi ()27)1 are zero. For 
instance 



(D 2 A)W ~ D 2xy 



(A y F u ) [(A u A u )-(B u B u )]+(A y A u ) [2(A U F U )+2(B U G U ) +Tr (xuXu 



= 
(28) 



and similarly 

~ (F X F U ) [(A U A U ) - (B U B U )] + (F X A U ) [2(A U F U ) + 2(B U G U ) + Tr(XuXu)] = . (29) 

The Feynman diagrams corresponding to the different contributions in ()29)1 are depicted in 
fig- 1. ' 

i \ i \ 

WvAAA) I 'WWW I 

\ I \ I 

\ / \ / 




/ W\, ? S WV f 1 'VW- 




Figure 1: Tadpole cancellation. The (bold) curly and (bold) dashed lines denote the auxiliary 
field (G) F and the scalar field (B) A, respectively; the solid line denotes the fermion field. 

The vanishing of the A and F one-point functions is due to the exact cancellation of 
the tadpole diagrams on the lattice. Similarly, the G and B one-point functions are zero 
at this order due to the presence of a matrix 75 inserted in the fermion loop. In order to 
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prove the WTi (|27|) one has to show that also the contribution depending on R vanishes. 
Indeed one finds 



(i?«>W = ((l-^ 1 )- 1 A + m)- 1 (AL J ,>(°) 



{XxXy) 



2{A y D 2yz A z ) - 2(B y D 2yz B z ) - D 2yz (A z A z ) + D 2yz (B z B z ) 



0, (30) 



where (fT3|) has been used. 



4 Two-point Ward-Takahashi identity 



In this section we discuss a more interesting WTi that relates the fermion and scalar two- 
point functions. Taking the derivative of with respect to fj and J a and setting to zero 
all the sources one obtains 



(XyXx) - (D 2yz (A z - i^B z )A x ) -{{F y - i^G y )A x ) + g{R y A x ) = . 



(31) 



Making use of the propagators given in (jl5j) . this identity is trivially satisfied at tree level. 

The next non-trivial order is g 2 which corresponds to the one-loop diagrams and can 
be written as 

{ Xv X x ) {2) -{D2y Z {A z -i lb B z )A x )^-{{F y -^ 

(32) 

where we used the expansion (JTBJ) for the function R. 

Applying the Wick expansion, the first term of this WTi is 

(XyX X } (2) = jiXyXx E [*( A + t^B)* + F ^ ~ + 2GAB 



X 



X (A + i l5 B)x + F(A 2 -B 2 ) + 2GAB )(°> . (33) 
We first isolate among the various contributions the tadpole ones 

(XyXx)? = g 2 Yl { fax*) (x*Xx) [(A Z F U ) {{A U A U ) - (B u Bj 

zu 

+ 2{A Z A U )[{A U F U ) + {B U GS) ~ (A z A u )Ti( X uXu) 

+ (XyXz)l5(XzXx)(B z B u )Tr(xulr J Xu)} ■ (34) 

Using the propagators (JT3J) and the relations Tr(x7 5 x) = and Tr(xx) — 4(AF) = 4{GB), 
it is easy to demonstrate that the tadpole contributions cancel out (see Appendix B). This 
property is general and also holds for the other terms of the WTi (I32jl . Therefore, one is 
left with the connected non tadpoles diagrams 



= 2g 2 ^2[(XyXz)(XzXu)(XuXx)(A z A u ) - (XyXz)l5(XzXu)l5(XuXx)(B z B u )j. 

uz 

(35) 

The corresponding Feynman diagrams are given in fig. 2. 
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-H5 



Figure 2: Non-tadpole contributions to (xx)^- 



The non-tadpole contributions to the second term of (|32j) are (here and in the following 
the sum over repeated indices z, u, w is understood) 

(D 2yz (A z - i^B z )A x ) { ^ T = g 2 \D 2yz {A z A u ) [Tr {{XuXw) (XwXu)) + 2{A U A W ){F U F W ) 
+ 2{B U B W ){G U G W ) + 2(F U A W )(A U F W ) + 2(B U G W ){G U B W )} (A W A X ) 
+ D 2yz (A z F u ) [(A U A W )(A U A W ) + (B U B W ) (B u B w )j (F W A X ) 
+ 2D 2yz (A z F u ) [(A U A W )(A U F W ) - (B u B w )(B u G w )j (A W A X ) 
+ 2D 2yz (A z A u ) [(A U A W )(F U A W ) - (B U B W ) (G u B w )j (F W A X ) j . (36) 

The corresponding Feynman diagrams are given in fig. 3. 







'VWj 



rVW 



-a/wi 



wv 




(WW ■ 



X" "s. 





Figure 3: Non-tadpole contributions to {D 2 (A — i^§B)A)( 2 \ 
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The non-tadpole contributions to the third term of (|32|) are 

({F y - ij 5 G y )A x ) {2) T = g 2 [2{F y A u ) \^([ Xu x w ){x w Xu)) + (F U F W )(A U A W ) 

+ (F U A W )(A U F W ) + (G U G W )(B U B W ) + (B U G W )(G U B W ) (A W A X ) 
+ (F y F u ) [(A U A W )(A U A W ) + (B u B w )(B u B w )j (F W A X ) 
+ 2(F y A u ) [(F U A W )(A U A W ) - {G U B W ){B U B W )\ (F W A X ) 
+ 2{F y F u ) [(A U A W )(A U F W ) - {B U B W ){B U G W )\ (A W A X ) 
- l5{G y B w )Tr(j 5 {x w Xu){XuXw)) (A U A X )\ (37) 
and the corresponding Feynman diagrams are given in fig. 4. 



/ WV 




wv-- 




/ vw- 




'VW- ■ 



O 



'WW 



rvw- 



VW\J 



'vw- 




)VWV - 




)rvw- 



/ w\/v( 




W\A| 





Figure 4: Non-tadpole contributions to ((F — i^G)A)( 2 ). 

Notice that the terms in last two rows of (JHHJ) cancel out since (AA) = (BB) and 
(AF) = (BG). These terms, originating from the last four diagrams in fig. 3, are the one- 
loop contribution to the 1PI AF-vertex function, which therefore vanishes at this order. 
Similarly, the last three rows of (|37[) do not contribute. In particular the last term, i.e. the 
last diagram of fig. 4 vanishes and that gives (G y A x )^ = 0. 

For the terms of the WTi (J32|) involving the function R one finds 

<iif = -J>{ X x} yz (AL z [ X (A + t l5 B) x + F(A 2 - B 2 ) + 2GAB]A X )^ , (38) 
where the fermion propagator follows from (|19|). 
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Also in this case the tadpole diagrams cancel out and one is left with 



x {2 



(A z F w )D 2zu (A u A w ) + (A z A w )D 2zu (A u F w ) - D 2zu (A u F w ) (A U A W ) 
— (B z G w )D 2zu (B u B. w ) — (B z B w )D 2zu (B u G w ) + D 2zu (B u B w ) (B U G W ) (A W A X ) 



+ \2{A z A w )D 2zu {A u A w ) - D 2zu (A u A w )(A u A w ) 

+ 2(B z B w )D 2zu (B u B w ) — D 2zu (B u B w ) (B U B W ) (F W A X ) \ . 
The corresponding Feynman diagrams are given in fig. 5. 



(39) 





rvw 



tm 



Figure 5: Non-tadpole contributions to {R^ A)^\ The blob denotes the insertion of the operator 
D 2 acting on the three legs outgoing from the vertex as in equation (|39|) . 

Finally, for the last term of (J32|) one gets 

(RWA.)® = -V2(xx) yz ((A z + z l5 B z )(xx) zw AL w A x )W 

= -2{(x v Xz)(XzXw) (A z A w )D 2wu (A u A x ) + (A W A X ) D 2wu (A z A u ) 

~ D 2wu (A z A u ) (A U A X ) 
- (XyX,z)lh{XzXw)lh (B z B w )D 2wu (A u A x ) + (A w A x )D 2wu (B z B u ) 



— D 2wu (B z B u ) (A U A X ) 

and the corresponding Feynman diagrams are presented in fig. 6. 



(40) 



Figure 6: Non-tadpole contributions to (R^'A) ^ . The blob denotes the insertion of the operator 
D 2 acting on the three legs outgoing from the vertex as in equation (|4U|) . 
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4.1 Calculation in momentum space 

In order to verify the WTi (J3"2*Jl we find convenient to work in the momentum space repre- 
sentation. 

For the fermion two-point function, the sum of the two diagrams in (135)) gives 

( X {p)x{q)) {2) = V(2tt)V(p + q) (d 2 { P ) - m(l - ^(p))) J G^ip, k)D 2 {p + k) 



x (D 2 (p)-m(l--D 1 (p)) 



(41) 



where 



G(p, k) = [M(p)((l ~ |A(p))] 2 - [M(k + p)((l - ±D x {k +p))] (42) 

and Di(p), D 2 (p) and M.(p) are the Fourier transform of the operators given in (|T2"j) and 

(USD- 

Similarly, the terms in (|36|) and (|37|) in momentum space write 
(D 2 (p)(A(p) - i lh B{p))A{q))^ = (D 2 (p)A(p)A(q))W = g 2 {2nf5\p + q)D 2 {p) 

x / g-\p,k)\2m\l - ^i(p)) 2 -Tr[D 2 (k)D 2 (p + k)] + 4D 2 2 (k)] (43) 



and 

((F(p) - i lhG {p))A(q))W = (F(p)A(q))W = mg 2 {2^5\p + q)(l - |a(p)) 
x JG~\V, k) [Tr(D 2 {k)D 2 (p + fc)J - 4D 2 2 (£ ; ) - 2 J D 2 (p) 
respectively. 

Finally, the two terms in (|39j) and (|4Uj) involving R in momentum space become 
(RP(p)A(q))V> = 2mg 2 (2K) 4 6 4 (p + q)(l - -D^p)) (d 2 ( P ) - m(l - ±D x {p)) 



(44) 



x jg-\p,k)( [ 2D 2 {p + k)-D 2 {p) 



(45) 



and 



(R {2 \p)A(q))^ = -Ag 2 {2^5\p + q) ( D 2 (p) - m(l - -D^p)) 



x J g~ l {p, k)D 2 (p + k) \D 2 {k) + D 2 (p) - D 2 (p + k) 



(46) 



In order to verify that the WTi (j32j) is exactly satisfied, we find convenient to arrange 
the various terms according to the powers of m. 

Inserting (gTJ) and flUD-flM)) into the WTi (jH21) and setting m = one has 



4g 2 (27r) 4 5\p + q) Q~\p, k) [D 2 (p)D 2 (p + k)D 2 (p) + D 2 (p) [D 2 {k) ■ D 2 (p + k) - D 2 2 (k) 



D 2 (p)D 2 (p + k) D 2 (k) + D 2 (p) - D 2 (p + k) 



(47) 
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where ^(7^7^) = 45^ has been used. Taking advantage of the invariance of Q(p, k) under 
the change of variables k — > —k—p, one can replace D 2 (p+k)D 2 (k) with | {D 2 {p + k), D 2 {k)} 
D 2 (p + k) ■ D 2 (k) and therefore the integrand exactly vanishes. 
The terms proportional to m add up to 



g 2 (2n) 4 5 4 (p + q)(l- -D x {p)) J Q~\p, k) [ - A[D 2 {p)D 2 {p + k) + D 2 (p + k)D 2 (p) 
+ AD 2 2 (k) + 2D 2 2 (p) - AD 2 (k) ■ D 2 (p + k) + £> 2 (p) (4D 2 (p + k) - 2D 2 (p) 



+ 4D 2 (p + k)[D 2 {k) + D 2 (p) - D 2 {p + fc) 



(4? 



Performing the substitution D 2 (p + k)D 2 (k) — > D 2 (p + fc) • D 2 (k) as described above it is 
easy to check that (|48jl vanishes. 

Finally, the contribution left is the one proportional to m 2 , i.e. 



g 2 W5\p + q)(l- -D^p)) 2 J Q-\p, k) \AD 2 {p + k) - 2D 2 (p) 

- 2{2D 2 {p + k) - D 2 {p)) 



(49) 



which is trivially zero. 

This end up our proof that the WTi (p?2*|) is exactly satisfied at finite lattice spacing. 



5 Continuum limit 



In this section we study the continuum limit of the WTi (JHzj) and discuss the restoration of 
the continuum supersymmetry in this limit. This will clarify the mechanism of cancellation 
between the different terms in the WTi and the role of the operator g(R(p)A(q))( 2 \ 
Following the notation of Ref. the operator D in © can be written as 



where 



and 



D(p) 



u(p) = - 
a 



j Eg 7m sin(jya) a 
2[u(p) + b(p)} + 2 



^sin 2 (p M a) + (ab(p)f 



1 1/2 



b(p) = ~ 



£2sin 2 (^)-l 



With this notation, the operators D\ and D 2 (D = D\ + D 2 ) are 

u(p) + b[p) 



Di(p) 



and 



D 2 (p) 



a?uj{p) 



auj(p) 



(50) 



(51) 



(52) 



(53) 



(54) 
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Similarly, 



and 



uj — b 
2u 



2iua 2 



(55) 
(56) 



Each term in the WTi 



as 



4(uj + b) + a 2 m 2 (uj - b) 

is a function of the external momenta p and can be written 
d A k 



Hp) 



(M 



rF(k, P ) 



(57) 



where the integration momenta k G [— -, -]■ If the integral (J57)) is ultraviolet convergent, its 
continuum limit is obtained substituting the function F(k,p) with its continuum equivalent. 
Otherwise, if (|BT|) is divergent and contains only massive propagators so that F(k,p) is finite 
for any set of exceptional momenta, one can use the lattice version of the BPHZ technique 
[TH] by writing 

Hp) = I C (P) + I\v) , (58) 



where 



f d^k \ 1 

np) = J (2^ [ F{klP) ~2Z-\P^---p^ 



Ml "Mn / p=0. 



dp, 



I\p) 



d A k 



n F 

Y- 

n=0 



P[tl ■ ■ -P/J-n 











dp, 



dp» n 



F(k,p) 



(59) 



(60) 



p=0 



and np is the degree of divergence of the diagram. I c (p) is ultraviolet finite and therefore 
its continuum limit can be taken. All the effects of the lattice regularization remain in 
I l (p), which is simply a polynomial in the external momenta with coefficients given by 
zero-momentum lattice integrals. In the following, we compute the lattice contributions of 
the Green functions entering in the WTi (J32*|) . Before doing this computation, we comment 
on their continuum part, such as I c (p), containing the subtracted integrand. Since the 
subtraction makes the integrals UV finite, the order of the limit of zero lattice spacing 
and the momentum integral can be interchanged. Applying this procedure to (RA) one 
immediately recognizes that its continuum part vanishes, since the function R vanishes for 
a —> 0. This is also clear due to the presence in and pp^) of the factor 1 D 2 (k) + 
Dzip) — F) 2 {p + k) that vanishes in this limit. For this reason one can restrict the analysis 
of the WTi to their lattice part. 

For the fermion two point function (|41|). one has to consider the following integral 



d A k 



1 



UJ 



p+k 



a J (2vr) 4 [(a/ + b') k + ^f{u' - b%] [{u> + b% +k + &f{w> - b% +k ] 
where k has been rescaled to k — > k/a and we have defined 



^7 M sin(£v) 



u' k = au(k/a) 



i-4E sin4 (y) + 4 (E sin2 (y)) 



2-i 1/2 



(61) 



(62) 



1 Actually, in (|45|l one must first make the change of variables k 
D 2 {p + k)- D 2 {k). 



-k — p to rewrite 2L>2(p + k) as 
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and 

uj' p+k = auj(p + k/a) 



1 - 4 £ sin 4 ( (fc + a ^ ) + 4 ( J] ^ ( * + ° P) ' 



2-i 1/2 



• (63) 



Similarly, o' fe = ab(k/a) and o' +fc = ao(p + k/a) and their espressions can be easily read 
from (JE2J)- 

The factor 1/a in (jfilj) implies a linear UV divergence of this integral which is cured by 
performing a Taylor expansion in pa up to the first derivative. The first term of the Taylor 
expansion of (J6*T]) is odd in k, thus is zero, while the first derivative is 



d 4 k ~* E M 7 M sin(/c M ) v-^ d 



with 



(27T)4[(^ + fo /) fc + £ 







dppd 



UJ 



p+k 



(64) 



dp p a 



UJ 



p+k 



[(uj' + b% +k + ^(oj'-b% +k } 



=0 [{uj> + b% + ^(uj>-b%} 



(65) 



plus terms proportional to a 2 m 2 which do not contribute in the limit a — > 0. Notice that 
in the denominator the term proportional to a 2 m 2 must be kept in order to ensure the 
IR finiteness of the integral, since (u/ + b')k ~ k 2 /2 for fc — > 0. Indeed, by substituing 
this derivative in (JMj) one sees that the contribution from the last term of ([65)1 produces a 
log(a 2 m 2 ) divergence (for a — > 0) originating from the k — > integration region, while the 
remaining terms give rise to a finite integral. 

Therefore, including the external leg factors, the fermion two point function can be 



written as 



where 



(xx) [2 \p) 



-C 2 i p- 



m) 



c 2 = i 



(p 2 + m?) (p 2 + m?) 

d 4 k to' k sin 2 (k p ) 



V 



(2K)*[( u + b y k + £f.(u>-by h ]* 

and C*2/ is a finite constant that, for our purposes, need not to be computed. 

For the scalar two point function (|4*3J) one has to calculate the following integral 



(66) 
(67) 



d A k { m 



UJ 



p+k 



(27r)H 2 [(uj' + b') k + ^(uj'-b%}[(uj' + b% +k + ^(uj'-b') p+k } 



UJ 



P+k 



a 2 {uj> + b% + ^f{uj' - b% (cu> + b') p+k + ^f{u> - b% +k 



X 



£,i sin 2 (/c At ) ^ sin(/c M ) sin(A; M + ap p ] 



<4 



p+fc 



(6$ 



The first term can be evaluated directly at pa = while for the second we need a Taylor 
expansion up to the second derivative in pa due to the factor 1/a 2 . We first concentrate 

2 From now on the factor (2tt) 4 S^ 4 ' (p + q) will be understood. 
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on the latter term. It vanishes at pa = and moreover its first derivative is odd in k and 
therefore also this term of the expansion vanishes. Thus the scalar two point function is 
given in terms of the following integral 



d A k f m 



[to 



I \2 



(2tt) 4 1 2 [^ + h ') k + ^f(u'-b') k ] 



( u ' + b') k + ^(cu'-b') k 2 



d 2 



pa 



dp p adp a a 



p+k 



p+k 



X 



Eu sin 2 (^) V sin(^) sin(^ + ap p ] 



p+k 



J p=0 



(69) 



There are two contributions coming from the second derivative. One is given by the product 
of f)I)J) with 



dp a a 



p+k 



sin(/c CT ) cos(/c CT ) 



2E M sin 2 (^)E^ CT sin 2 (^)sin(A ;CT 



(CO 



M3 



(70) 



which produces a log(a 2 m 2 ) divergence to (JBTJj) . originating from the product of the last 
term of (jB3j) with the first term of (fTUj) . The second is given by the second derivative of 
the third line of (J69|) and its explicit expression is not needed since its contribution to the 
integral (J69j) is finite for a — > 0. 

Collecting all terms and including the external leg factors, the two point function 
becomes 



D 2 {AA) { - 2 \p)=i £, 



(p 2 + m 2 ) (p 2 + m? 



-Cm 2 - C lP 2 



where 



C 3 = g> 



d 4 k 



{2-KY[^ + b') k + ^f(uj'-b') k } 2 
d A k sin 2 (k p ) cos(kp) 

(27T) i [(u' + b') k + ^(uj'~b%f + 



c 



If 



(71) 

(72) 
(73) 



and Cif is a finite constant. 

A similar analysis applied to (|44|) and gives 



(FA) (2) {p) = m 



(p 2 + m 2 ) (p 2 + m 2 ) 2 



C 3 + d)p 2 



(74) 



The continuum limit of the two point function containing the operator R can also be 
determined. For (|45|) and (^HJ) one has 



(R w A) w (p) =m 



m) 



-5— — 5t(C2 — -^C?)i 
(jr + m z ) (p z + m" 1 ) 2 



and 



(rWa)®(p) 



(i fa — m) 



(p 2 + m 2 ) (p 2 + m 2 
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■(^-dy 



(75) 



(76) 



respectively. Notice that the combinations C9 — C\ and C2 — \C 3 are two (different) finite 
numbers. Indeed, looking at the k — > behavior of the integrand of (|67|). (|72|) and (J75j) 
one sees that the log(a 2 m 2 ) contributions cancels out in these combinations. This is a 
consequence of the fact that the one-loop correction to the two-point functions of A, F and 
X have the same logarithmic divergent parts [20J. 

Substituting Ijfifij) . (JHJ) and ((7i)l - (|7K)l in ({HZ)) one verifies this WTi in the continuum 
limit : 

(1 



(p 2 + m 2 ) (p 2 + m 2 ) 

"U m C 3 —p Ci)- 



[p 2 + m 2 ) 2 (p 2 + m 2 ) 

■(Ci + -C 3 )p 2 



(p 2 + m 2 ) 2 (p 2 + m 

(i 6 — m) ,, . w _ 1 _ 



2\ 



(p 2 + m 2 ) 2 (p 2 + m 2 ) 

+ i ff~ 7 t ( C2 ~ C ^ 2 TT^S = • (77) 
(jr + rrr) (jr + 777/J 

Actually, this is a check of the results we have obtained for the continuum limit of the two 
point functions, since we have proved that this WTi holds for any a and therefore must be 
verified also in the limit a — > 0. Notice that the term (RA) in (|31|) is essential to recover 
the WTi 022) also for a -> 0. 

Let us clarify the role of the operator R. Thanks to the exactness of WTi (|H2*|) it is 
always possible to write the two point function (RA)^ as a suitable combination of the 
other three two point functions involved in this WTi. In particular, in the continuum limit 
one can write 

\ RA ) = 9 % P^^T, 9 + 1 9 (^P + o 3 m 



pZ _|_ m z pZ _|_ m z pZ _|_ m z pZ _|_ m / 

Vil , „ „ , n 1 



(*a - 5 3 )p'^—- (71 



p2 _|_ m 2 ^2 _|_ m 2 

where, from (J7HJ) and (f7Hjl . 

5i = ic 3 -C 2 -5 3 , 52 = ^3-^-53, (79) 

and the constant £3 is arbitrary. Then in the continumm limit one can rewrite the WTi 
(jH2} as the supersymmetric continuum WTi 

(XX) { S - i i>{AA)% - (FA)f = (80) 

with 

m , „ i i> — m 



(xx)g } = (xx) (2) + ^P?^ A 

+ mr p z + m z 
(AA)%> = (AA)(2) - -^—^P 2 + hm 2 )^— 2 

\ ) pZ _|_ m z pZ _|_ m z 

{ FA)$ = (FA)V + -^(8 2 - 8 3 )p> 



pi _|_ v - ^2 _|_ m 2 
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It is convenient to express these two point functions in terms of 1PI vertex functions: 



i to — m ro\ i to — m 



p 2 + m 2 xx p 2 + m 2 



pZ _|_ m z pZ _|_ m z 

(FA)^ = -^T-^aa ~ P 2 Z%)^~ 2 , (82) 
p l + m z p l + m 2 

where the vanishing of the 1PI AF- vertex has been used. ^From (|66|). (|7T|) and ()74|). the 
lattice contribution to these 1PI vertices in the continuum limit reads 

Eg = i#7 a , E^^^d, Sg>. = -|ck . (83) 

Moreover, from (|81|1. one has 



= i fa 


c 2 , 


T m _ 


= p 2 c, , 


! has 








.(2) 
'XX B 


- V (2) 


f i j6<5i 




.(2) 
'AAA 


- V (2) 


+ p 2 5 2 


= p( t - 


.(2) 
'FFR 


- V (2) 
— FF 


+ 5 3 = 


-(— - * 
v 2 



-^AP 2 



54) 



with 



Z x = Z A = Z F = -(^-5 3 ). (85) 

In Ref. j2D] it was shown that the one-loop corrections to the two-point function of A, 
F and x differ by finite quantities. Our construction shows that if one redefines the 1PI 
vertices as in f)84j) the wave function renormalization factors become equal. This is an 
important consequence of the exact lattice supersymmetry we have introduced and of the 
WTi derived from this symmetry. This automatically leads to restoration of supersymmetry 
in the continuum limit with equal renormalization wave function for the scalar and fermion 
fields. 

In a more standard approach [T^ I2T] the function R is not included in the lattice 
supersymmetry transformation. Since the action is not invariant under this transformation, 
the WTi contains a breaking term. From the a — > limit of this WTi one determines the 
counterterms nedeed to restore supersymmetry in the continuum limit. The central issue of 
our approch is that this possibility is guaranteed by the existence of an exact supersymmetry 
of the lattice action. 



6 Conclusions 

In this paper, starting from the N = 1 four dimensional lattice Wess-Zumino model that 
uses the Ginsparg- Wilson relation and keeps an exact supersymmetry on the lattice, it is 
showed that the corresponding Ward-Takahashi identity is satisfied, both at fixed lattice 
spacing and in the continuum limit. This result crucially depends on the Ginsparg- Wilson 
properties of the operators involved in the lattice action. The calculation is performed in 
lattice perturbation theory up to order g 2 in the coupling constant. 

It is also showed that the study of the continuum limit of this Ward-Takahashi identity 
determines the finite part of the scalar and fermion renormalization wave functions which 
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automatically leads to restoration of supersymmetry in the continuum limit. In particular, 
these wave functions coincide in this limit. 

Although we limit our computation up to the order g 2 , this order is not trivial and the 
discussion is general and applies to higher orders by following the procedure described in 
this work. 

There are several issues that remain to be investigated. First of all it will be interesting 
to perform numerical simulations of this model to check non-perturbatively the WTi (J31I) . 
Furthermore, one of the most important question is whether these ideas may be extended 
to theories with a gauge symmetry. 
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7 Appendix A 

The fermionic kinetic term in the action (j2J), 

(1 - -D^D, (86) 

needs a careful study at the border of the Brillouin zone. To this end we set 

p„= (0,0,0,—) (87) 
a 

and study the limit e — > 0. Using Eqs. (j51|) and (J52j) we have that 

b = i[l - t + {e")] (88) 

and 

lo = - . (89) 
a 

Inserting these values in and we find 

(l-^)-^- [1 + 0(6*)] (90) 

and 

D 2 = i l4 -[l + 0(e 2 )]. (91) 
Finally, the fermionic kinetic operator behaves as 

(1 - %D X )~ X D 2 = [I + 0{e 2 )] (92) 

thus in the limit e — ► with fixed a the would be doubler (|57jl becomes (infinitely) massive. 
Similarly, one can check that this value of the momentum do not generate a pole in the 
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bosonic propagators. This analysis can be generalized to the other edges of the Brillouin 
zone. For instance, if 

7r — e 7T — e. . 
p„ = (0,0, , ) (93) 



we have 



CL 3 

{l-^-D 1 )- 1 D 2 = i( l3 + lA )- 



(94) 



which again becomes a massive mode when e — > while a is kept fixed. It is easy to 
demostrate that all the rest of the would be zero modes behave in the same way. 
Notice that the fermion propagator in (|15|) can be rewritten as 



(XX.) = ~(D 2 - m(l - -Dx)) 



-Di + m (1 - -Di) 

CL — 



(95) 



which is clearly finite for e — ► 0. 

Appendix B 

In this appendix we explicitly show that the tadpole contributions to the two point WTi 
(132)) cancel separately. For the (xyXx) two point function this has been already shown 
in Section 5. 

The tadpole contribution to (D 2yz (A z — i'~f 5 B z )A x ) is 

(D 2yz (A z - i^B z )A x )f = g 2 {D 2yz (A z A u ) [(F U A W ) (2(F W A W ) + 2{G W B W ) - Tr( Xw X w ) 

+ (F U F W )((A W A W ) - (AA))] (A U A X ) 
+ D 2yz (A z F u ) [(A U A W ) (2{F W A W ) + 2(G W B W ) - Tr{ Xw Xw) 
+ (A U F W ) 

+ D 2yz (A z A u ) [(A U A W ) (2(F W A W ) + 2(G W B W ) - Tr( Xw X w ) 

+ (A U F W ) ({A W A W ) - (B W B W ))] {F U A X )Y (96) 

Similarly, the tadpole contribution to ((F y — i^G y )A x ) is 

((F y - i^G y )A x )f = (F y A u ) [(F U A W ) (2(F W A W ) + 2(G W B W ) - Tr( Xw Xw)) 

+ (A U F W ) y(A w A w ) — (B W B W )\ (A U A X ) 
+ (F y F u )(A u A w ) [2(F W A W ) + 2(G W B W ) - Tr( Xw Xw)] (A U A X ) . (97) 
It is easy to see that both expressions are exactly zero. 
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